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Abstract 

In a realistic application of the SPA + RPA theory for calculation of the 
nuclear level densities we find that quantal fluctuation corrections (RPA) 

are important even up to temperature T = 2.0 MeV. This leads to a good 

• i— i . 

/\ ' agreement between calculated numbers and the available experimental data 

H : 

for 104 Pd and 114 Sn, particularly the excitation energy (E*) dependence. 
Furthermore, we also argue that a = S 2 /4E* is the only correct definition 
of the level density parameter in the present context which is also consistent 
with the Bethe like level density formula. 
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In the last a few years there have been considerable efforts to develop micro- 
scopic methods for the calculation of accurate values of level densities as a function 
of excitation energies]!]]- 1|. One of these methods which is of present interest is 
based on the auxiliary field path integral representation of the partition function for 
a given nuclear Hamiltonian. The path integral representation can be obtained in 
two ways: (a) the so called as shell model Monte Carlo (SMMC) method[(|, H and 
(b) the SPA+RPA approach which includes the thermal fluctuations through 
static path approximation (SPA) and the quantal fluctuations about static paths 
are included using random phase approximaton (RPA). It is now well established 
that SMMC approach can be applied even at very low temperatures. On the 
otherhand, SPA+RPA approach is computationally faster than SMMC approach 
and can be used for moderately low to high temperatures, typically, for T > 0.2 
MeV. However, in most of these studies so far the main emphasis has been to 
demonstrate, through simple model studies, the relative defferences between the 
SPA and SPA+RPA results for the temperature dependence of the energy and the 
level density. 

Recently f7|], the level densities as well as the level density parameters have been 
extracted for the medium mass nuclei 104 Pd through the measurement of proton 
yields in reaction 93 Nb( 12 C,p) W4 Pd and that of lu Sn through 103 Rh( 12 C,p) lu Sn 
reaction. These data are available for the excitation energy (E*) ranging from 5—25 
MeV or equivalently T ~ 0.5 — 1.5 MeV which is well suited to test the feasibility 
of the SPA+RPA approach. In this letter we calculate the level density and the 
level density parameter as a function of E* for 104 Pd and nA Sn using SPA+RPA 
approach with a quadrupole-quadrupole interaction model Hamiltonian 

H = H - \ X E (Q,) 2 . (1) 

A fi=-2 

In the above, H represents the spherical part, Qq = Q' , Q +fl = -m{Q'^ + Q'£) and 
Q-v- = TF^Q'v ~ Q/t) w ith /i = 1 and 2 and Q's stand for the usual quadrupole 
moment operators. Value of the quadrupole interaction strength \ = 120A _5 / 3 / c 



MeV (A denotes the mass number) is taken from Ref. M where f c = 1 — 2 is a 
core polarization factor. 

The grand canonical partition function in SPA+RPA takes the following form 



Zrpa = 4vr 2 (^) 5/2 / P 4 d[3 J | sir&y | dje'^ 
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where, a = (huo) 2 /x with hcu = A1A x / 3 MeV . The quantities H' and Crpa are 
the single-particle Hamiltonian and the RPA correction factor, respectively, given 
by 

H' — H - huJofl (Q cos7 + Q +2 sin^) (3) 



and 



with 



Crpa =-- ( II Det I ° r ' 
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where, /y = f% — fj and Ay = 6j — e,- with /, being the Fermi distribution function 
and €i is the eigenvalue of H'. In the above, | i) represents an eigenstate of H' . 
The grand canonical trace in eq. @) can simply be performed using 



lit 1 + e- p ' e ' +a "]\ (nil + e -P'^ +a "}\ 



Tre~ ffH = lit 1 + e P \ IK 1 + e J ] ( 6 ) 

where, (3' = 1/T and a p (a n ) is the Lagrange multiplier required to adjust the 
proton (neutron) numbers. 

The SPA representation of the partition function can be obtained by putting 
Crpa — 1- It is, therefore, clear from eqs. (|4],|) that for higher temperatures 
Crpa — ► 1 or in other words Z RP a — ► Z S pa- 

Once the partition function is known, the level densiy or more precisely the 
state density W(E) can be calculated using saddle point approximation which 



gives 



F S 



where, 

S = InZ + f3'E - a p N p - a n N n (8) 

is the entropy with Z being the grand canonical partition function in the SPA or 
SPA+RPA approach. The quantities a P)U and (3' (or T _1 ) are so chosen that the 
saddle point conditions 

E = -p nZ (9) 

and 

N p , n = -= InZ (10) 

oa P:n 

are satisfied. The quantity T> in eq. (7) is a determinant of 3 x 3 matrix defined 
by the elements 

d - = di^ lnZ (11) 

where, x = (/3', — a p , —a n ) and the second derivatives are evaluated at saddle 
points. For the sake of comparison with the experiment we present below the 
results for p(E), instead of W(E), which can be obtained as[p"0|. 



where, a 2 = I ng jU is the spin cut-off factor with I rig being the rigid - body value 
of the moment of inertia. 

The model space used to perform the numerical calculations is given in Table 
1. For the range of temperature T < 2.0 MeV this basis space should be adequate. 
The choice of the values of the spherical single particle energies is a rather difficult 
task. Due to the quantal nature of the nucleus there is no smooth A-dependence 
for a large range of A. From our experience in the p/-shell and rare earth region we 
have finally chosen these numbers with the help of Fig. 1 (suitable for A p» 100 ) in 
Ref. |11[ and Fig. 1 (suitable for 108 Sn) in Ref. |12|] which are actually obtained as 



solutions of Woods-Saxon potential. The neutron core with N=40 ensures sufficient 
number of active valance neutrons. The value of the core polarization factor, f c , is 
taken to be 1.5 . With this the value of the quadrupole deformation parameter, j3, 
in the ground state of l0A Pd comes out to be about 0.1 which is quite reasonable 
[ PI . lu Sn is spherical in the ground state. The pairing correlations are also 
expected to be small for these nuclei. 

For the sake of compactness most of our results and discussions will be presented 
for 10A Pd only. For 11A Sn results on level densities will be presented towards the 
end , before conclusions. In Fig. 1 we display the SPA as well as SPA+RPA 
results for the variation of energy as a function of temperature for l0A Pd . To 
obtain E(T) within the SPA+RPA we first of all check the convergence of the 
RPA correction factor C- 1ZVA by choosing various values of N m in eq. (4). We find 
that N m = 40 is sufficient but we have used N m = 80 in the present calculation so 
that at very low temperatures it should be sufficiently accurate. We see that, as in 
Refs. H [14|], the RPA or quantal fluctuation corrections lower the value of E(T) 
at lower temperatures. As temperature increases, quantal fluctuation decreases 
or C nvA — > oo yielding the value of E(T) close to the one obtained within SPA. 
However, it shows that the RPA corrections are important up to about T = 2 
MeV. We then next show in Fig. 2 the variation of the level density as a function 
of E*. The values of E(0) needed to calculate E* are obtained in the case of 
SPA as well as SPA+RPA by extrapolating the corresponding curves of Fig. 1 
to T = 0. These values come out to be equal to —310.603 and —315.357 MeV 
in the case of SPA and SPA+RPA, respectively. We see that SPA+RPA results 
for the level density are in good agreement with the ones extracted recently [|7| 
through the measurements of proton yields in the reaction 93 Nb( 12 C , p) 104 Pd. On 
the other hand, SPA level densities are higher compared to the measured values. 
For instance, the ratio Pspa/ pexp ~ 10 2 at E* = 10 MeV which reduces to about 
10 at E* = 24 MeV. Whereas, Pspa+rpa/ 'pexp varies from 1.0 — 1.2 over the entire 



range of E* for which experimental data are available. It should be mentioned 
that we have not put the error bars on the experimental data. The experimental 
data have errors of the order of 20 - 30% at the lower end and of a few percent at 
the other end. 

Next, we now discuss about the parameterization of the SPA+RPA level densi- 
ties in order to extract a value of the level density parameter, a. We have shown in 
our earlier work[|15j that SPA level density can be reproduced by Bethe's formula 
provided an appropriate value of the parameter 'a' is used. Normally there are two 
relations which are used to compute the value of a, namely, 

E* = a e T 2 (13) 

and 

S = 2a s T (14) 

where, the suffices e and s are used to distinguish the value of 'a' determined using 
E* or S, respectively. Treating a e and a s equal, usually a third relation a = S 2 /AE* 
is also derived. However, in reality there is no rigorous reason to do so. Usually all 
the three relations yield different values of a when computed numerically. Hence, 
it becomes a natural question to ask as to which 'a' is appropriate for the use in 



Bethe's formula for the level densities 10 



P(E) = C'; , J^r (15) 
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It is important to realise from eqs. (7 - 9) that the constant part of the 
p ref actor fl~6|| in eq. (2) which normalizes the 'measure' does not contribute to 
the calculation of E* but certainly contributes to the value of S. In any mean field 
approach the effect of this prefactor would be missing. Now using eqs. fll3|) and 
( |I4"D we can write 



2/?* 



lalE 



S = 2J^ — = 2VaE* (16) 



where, 

a = a 2 Ja e = S 2 /AE* (17) 

We find that the values of a required in eq. QT5| ) to get p SP a+rpa is quite close to 
the one given by eq. fll7|) . For instance, at T = 0.5, 1.0, 1.5 and 2.0 MeV we 
get a(a fit ) = 10.46(11.04), 12.92(13.12), 12.42(12.49), 11.89(11.94), respectively. 
In Fig. 3 we have displayed the variation of the inverse level density parameter 



(K = A/a) as a function of the excitation energy per particle, e = E* /A [|17 | 
which correspond roughly to T = - 2.0 MeV. The values of K e , K s and K es are 
obtained using 'a' calculated from eqs. flT3]), (JT4J) and QTTD, respectively. We see 



that in comparison to K e and K s the values of K es are increasing very slowly with 
e or temperature for T > 1.0 MeV. Also, we would like to point out that the 
values of K es are quite close to the experimental values [0] available for E* = 16 
— 22 MeV. Behaviour of K in the very low temperature region(T < 0.5MeV) is 
reflecting the well known effects of shell structures. 

Finally in Fig. 4 we have displayed the variation of the level density as a 
function of E* for U4 Sn. The theoretical curves are drawn after dividing the 
actual numbers by a factor 20.827 such that the upper most point (with mini- 
mum uncertainty) of the experimental data matches with the one calculated in the 
SPA+RPA approach. The variation is rather very well reproduced by the solid 
curve keeping in mind the fact that the experimental values have also large inher- 
ent uncertainties J7|. The SPA curve is roughly parallel to the SPA+RPA one (the 
ratio Pspa/ Pspa+rpa ~ 4.5 at E* = 5.0 MeV and 2.6 at E* = 30 MeV) emptying 
that here the variation is approximately reproduced. 

In conclusion, we have studied the excitation energy dependence of level den- 
sities for l04 Pd and ll4 Sn including thermal as well as quantal fluctuations of the 
nuclear quadrupole shape parameters. We find that inclusion of quantal fluctu- 
ations is essential to reproduce the experimental data even upto the excitation 
energy of about 25 MeV (or T rj 1.5 MeV). We have also shown that the value of 



the level density parameter 'a' to be used in the Bethe's formula (eq. (|T5|)) should 
be computed from eq. (|17"D only. Also, this value of a can not be used to relate 
temperature to get S or E*, particularly in the SPA or SPA+RPA approach. Ac- 
tually the calculation of a from a 2 s /a e is essentially equivalent to fitting the value 
of the level densities to the Bethe's formula. 
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Table 1: Spherical 

N = 40 as a core for 



single-particle energies (in units of Huiq) with Z 
proton and neutron, respectively. 



28 and 



Protons 


Neutrons 


Spherical 


Energy 


Spherical 


Energy 


orbits 


(huj ) 


orbits 


(HuJo) 


1P3/2 


-1.376 


Ofi-9/2 


-0.975 


O/5/2 


-1.374 


14/2 


-0.484 


lPl/2 


-1.171 


Ofi-7/2 


-0.30 


0#9/2 


-0.975 


2si /2 


-0.216 


14/2 


-0.484 


14/2 


-0.122 


Ofi-7/2 


-0.30 


0/ill/2 


-0.122 


2S1/2 


-0.216 


04/2 


0.358 


14/2 


-0.122 


1/7/2 


0.405 


0/lll/2 


-0.122 









Figure Captions 

Figure 1: Variation of energy as a function of temperature. 

Figure 2: Variation of the level density as a function of excitation energy. Solid 
circles represent the experimental values taken from Ref. [7]. 

Figure 3: Variation of the inverse level density parameter (K = A/ a) as a function 
of e = E* / A. The curves labelled K e , K s and K es are obtained with the values of 
'a' calculated using eqs. (13), (14) and (17), respectively. 

Figure 4: Excitation energy dependence of level density for 114 S*n. Dashed and 
Solid curves represent the values of SPA and SPA+RPA level density, respectively 
reduced by a factor of 20.827 (see the text for details). 
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